Finite-size scaling in complex networks 
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A finite-size-scaling (FSS) theory is proposed for various models in complex networks. In par- 
ticular, we focus on the FSS exponent, which plays a crucial role in analyzing numerical data for 
finite-size systems. Based on the droplet-excitation (hyperscaling) argument, we conjecture the val- 
ues of the FSS exponents for the Ising model, the susceptible-infected-susceptible model, and the 
contact process, all of which are confirmed reasonably well in numerical simulations. 
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Critical phenomena in complex networks have at- 
tracted much attention recently and traditional models 
in statistical physics have been examined on diverse net- 
works 0, 0, y, |j, [E S • Interesting questions in such 
probes would be the existence of phase transitions and 
the network-dependent critical behavior near the tran- 
sition. Up to now, various equilibrium systems such as 
Ising, Potts, and XY models [J 0, EJGl as weu as nonequi- 
librium systems such as percolation 15 1, directed percola- 
tion and synchronization models [7| have been studied 
by means of the mean-field (MF) approach Q, the replica 
method 0], and the thermodynamic potential hypothe- 
sis It is predicted that phase transitions in com- 
plex networks exhibit mostly the standard MF critical 
behavior except in highly heterogeneous scale-free (SF) 
networks where the interesting heterogeneity-dependent 
(but still MF-type) critical behavior appears [l|, @, [1| • 

The MF prediction has been brought up to the test by 
extensive numerical simulations and passed it reasonably 
well in cases of less heterogeneous networks like random, 
small- world, and even some SF networks @,0,@,0|- How- 
ever, for the highly heterogeneous networks (like the most 
of SF networks found in nature) , the asymptotic scaling 
regime could not be reached easily in numerical tests due 
to huge finite size effects and consequently there is no rea- 
sonably solid numerical analysis reported as yet. There- 
fore it is essential to understand the finite-size-scaling 
(FSS) behavior analytically in networks, not only to an- 
alyze numerical data for finite-size networks but also to 
explore the physics of correlated size scales in networks. 

In the present Letter, we propose a FSS theory for 
various models in complex networks, based on a droplet- 
excitation (hyperscaling) argument. Our conjecture for 
the FSS exponent values is confirmed via numerical sim- 
ulations for the Ising model and the contact process Q . 

We first start with the standard FSS theory in low 
dimensional systems. As a typical example, consider the 
ferromagnetic Ising model. Its critical behavior near the 
transition is characterized by the singular behavior in the 
magnetization to ~ e' 3 , the susceptibility \ ~ l e l~ 7 > an d 



The standard FSS theory for the singular part of the 
free energy / reads 0] 

f(e, h, L^ 1 ) = b- d f{by T £ , V* h, bL~ x ), (1) 

where b is the scale factor, d the spatial dimension, L the 
system linear size, and h the external field. The two scal- 
ing dimensions, y T and y H , determine all thermodynamic 
exponents such as f3 = (d—y H ) /y T and 7 = (2y H —d)/y T . 
The FSS behavior of the spontaneous magnetization is 
derived from Eq.fT]) as 



= L- 0v Tg{L y T t ), 



(2) 



where g is the scaling function and d/y T =2/3 + 7. If 
the correlation length scales as £(c,h) = b£(b y T e,b v H h), 
the hyperscaling relation appears as y T — 1/v In 
this case, the FSS variable L y re can be expressed as 
(L/£) y r, implying that the correlation length competes 
with the system size. In low dimensions where the hyper- 
scaling holds, the above FSS theory works perfectly well 
for most equilibrium models and even for nonequilibrium 
ones (like directed percolation models) with appropriate 
modifications [Io| . 

In high dimensions where the MF theory becomes 
valid, the FSS becomes a little bit tricky. It has been 
argued that a "dangerous irrelevant variable" comes into 
play in high dimensions which breaks the hyperscaling 
(y T ^ 1/v) and modifies the FSS [HI]. Nevertheless, for 
block-shape samples with periodic boundary conditions, 
Eqs.fl]) and (F2J) are still valid with the trivial MF values 
of y T — d/2 and y H = 3d/4, which correspond to the 
MF bulk exponents such as j3 = 1/2 and 7 = 1. The 
broken hyperscaling implies that there must be another 
length scale competing with the system size, rather than 
the usual Gaussian correlation length. 

Consider the Landau free energy of the order parame- 
ter m near the transition in the MF regime as 



/(to) = —em 2 + um + C(m 6 ), 



(3) 



the correlation length £ ~ |e| 1 
temperature defined by e = (T c 



where e is the reduced 
- T)/T c . 



where the Gaussian spatial fluctuation term (Vto) 2 is ig- 
nored. In the ordered phase (e > 0), the free energy has 



2 



a minimum at m = m* ~ e 1 / 2 with / ~ — e 2 . One may 
ask a typical size £ T of a disordered droplet excitation 
out of the uniformly ordered environment. As the free 
energy cost by the droplet excitation is compensated by 
the thermal energy, (A/)£^ ~ fcgT, we find £ T ~ e~ u T 
with v T — 2/d — l/y T . The Gaussian length scale di- 
verges as £ G ~ e~ v o with u G = 1/2. 

For d > d u = 4 (the upper critical dimension), £ G 
dominates over £ T , which leads to the correlation length 
exponent v = v G and the MF theory is valid. However, 
the FSS variable L y re becomes (L/£ T ) y T, implying that 
the competing length scale is not the dominant correla- 
tion length but the droplet size. Substituting the linear 
size L by the volume N ~ L d , Eq. ((2]) reads 



m = Ar-<3/^(iV 1/p e), 



(4) 



where the FSS (droplet volume) exponent v = dv T = 2 
in the MF regime. For the general <p q MF theory (/ = 
—em 2 +um q ), we find that v = d u v G with d u = 2q/(q—2), 
which is consistent with the earlier result by Botet et al. 
for models with infinite-range interactions [l2j. 

We are now ready to explore the FSS in networks. Net- 
works have no space dimensionality and may be consid- 
ered as a limiting case of d — * oo. So we expect that any 
model in networks displays a MF-type critical behavior. 
In particular, the MF FSS exponent T> is independent of 
d, which leads to the natural conjecture that Eq.Q also 
applies in networks. These predictions have been con- 
firmed by numerical simulations for various models in 
random networks, small-world networks, and complete 
graphs. Moreover, the relation of v = d u v G has been 
exploited to calculate the value of d u via simulations in 
networks for complex nonequilibrium models [I3j . 

In SF networks with the degree distribution P(k) ~ 
fc~ A , there appears a nontrivial A-dependent MF critical 
scaling for A < X u (highly heterogeneous networks) while 
the standard MF theory applies for A > \ u 1]. Nat- 
urally, we expect a nontrivial FSS theory associated to 
the nontrivial MF scaling for A < X u . Previous studies 
pay attention to the MF analysis in the thermodynamic 
limit and hardly discuss the FSS in the general context. 
Recently a few numerical efforts have been attempted to 
confirm the MF predictions, but huge finite-size effects 
and the lack of the FSS theory disallowed any decisive 
conclusion for highly heterogeneous networks [J, [7J . Most 
recently, even a non-MF scaling has been claimed for the 
contact process [3 HH and a question arises whether the 
cutoff in degree k influences the FSS. 

We start with the phenomenological MF free energy 
for the SF networks proposed in [H, 12] 

f(m) = -em 2 + urn 4 + u|m| A-1 + C(m 6 ), (5) 

where the A-dependent term originates from the singular 
behavior of the higher moments of degree in SF networks. 
For A > A. u = 5, the A-dependent term is irrelevant and 



we recover the usual 4 MF theory, yielding (3 = 1/2 and 
v = 2. For 3 < A < 5, the A-dependent term becomes 
relevant and we find the (f> q MF theory with q = A — 1. 
A simple algebra leads to f3 = 1/(A — 3) and the free 
energy density in the ordered phase is / ~ — e 1+2/3 . One 
can estimate the typical droplet volume N T ~ (A/) -1 , 
yielding N T ~ <T» with v = 1 + 20 = (A - 1)/(A- 3). By 
including the external field term hm in Eq.([5|), one can 
show 7 = 1 for all A > 3. 

The results for /3, 7, and the FSS exponent v for the 
Ising model in SF networks are then summarized as 
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for 3 < A < 5, 
for A > 5. 



(6) 



For A < 3, no phase transitions occur at finite tempera- 
tures and, at A = 5, a multiplicative logarithmic correc- 
tion is expected [l[ . It is interesting to notice that a naive 
power counting for the <p q local theory with the Gaussian 
spatial fluctuation term (Vm) 2 yields the same result for 
v by using the relation of 9 = d u v G [ijj]. Our conjecture 
for v bears no reference to the degree cutoff k c caused 
by the finite system size N. We will argue later that the 
cutoff is irrelevant if it is not too strong: k c > N 1 ^ [151 ] . 

We check our conjecture via numerical simulations. 
Two typical SF networks are considered, namely the 
static model 1 171 and the uncorrelated configuration 
model (UCM) [181 ]. As these networks have different de- 
gree cutoffs (natural cutoff k c ~ N 1 /^- 1 ) versus forced 
sharp cutoff k c ~ iV 1 / 2 ) in finite systems, one may look 
for a possibility of the cutoff-dependent FSS behavior if 
any. It turns out that both cutoffs are not strong enough 
to influence the FSS for A > 2. 

We performed Monte Carlo simulations at various val- 
ues of A up to N = 10 7 . We measure the magnetization 
m, the fluctuation x' = -/V(Am) 2 , and the Binder cumu- 
lant B and average over ~ 10 3 network realizations. The 
transition temperature T c is estimated by the asymptotic 
limit of the crossing points of B for successive system 
sizes as well as of the peak points of \' ■ At critical- 
ly, Eq.® leads tom~ N~ / p and similarly x' ~ ^ 
with ~ |e|~ 7 in the thermodynamic limit. This power- 
law behavior in N provides an alternative check for the 
criticality as well as the estimates for the exponent ratios. 
In equilibrium systems, the fluctuation-dissipation theo- 
rem guarantees 7' = 7. By collapsing the data over the 
range of temperatures, we estimate the value of the FSS 
exponent D. Our numerical data for m and x' collapse 
very well for all values of A in both static and UCM net- 
works. In Fig. [Tl the data collapse is shown for A = 3.87 
in static networks. We summarize in Table I the numer- 
ical estimates for 0/v, and *f /v at various values of A 
in static and UCM networks. All data agree reasonably 
well with our predictions. 

We also measure the degree- dependent quantities like 
the magnetization on vertices of degree k, rrik, and its 
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FIG. 1: (Color online) Data collapse of m for the static 
network with A = 3.87, using j3/v = 0.37 and v = 3.5. Insets: 
double logarithmic plots of the critical decay of m and \' 
against N with slopes f3/u = 0.37(5) and j'/u = 0.26(4). 



TABLE I: Critical exponents of the Ising model on the static 
and the UCM networks, and the CP on the UCM networks, 
compared with our MF predictions. 







Ising 




network 


A 


p/u V 




MF 


A > 5 


1/4 2 


1/2 




3 < A < 5 


1 A-1 
A-1 A-3 


A-3 
A-1 


static 


7.08 


0.26(4) 2.0(2) 


0.45(5) 




4.45 


0.28(2) 2.4(2) 


0.45(3) 




3.87 


0.37(5) 3.5(3) 


0.26(4) 


UCM 


6.50 


0.24(4) 2.0(2) 


0.51(5) 




4.25 


0.31(1) 2.5(1) 


0.39(1) 




3.75 


0.38(6) 3.9(2) 


0.24(3) 






CP 




network 


A 


0/v V 


i/v 


MF 


A > 3 


1/2 2 







2 < A < 3 


1 A-1 
A-1 A-2 


A-3 
A-1 


UCM 


4.0 


0.49(1) 2.1(1) 


0.00(5) 




2.75 


0.58(1) 2.4(1) 


-0.16(2) 




2.25 


0.78(1) 4.0(5) 


-0.55(5) 



fluctuation (Amfc) 2 . These quantities are found to sat- 
isfy a scaling relation with the scale variable kN^ 1 ^ (not 
shown here) [la]. For k > N 1 '*, the distribution P(k) 
becomes almost flat for each realization of networks and 
the degree exponent A loses its identity. Therefore, ver- 
tices of such a high degree contribute in a trivial way 
and the cutoff beyond this range (k c > N 1 ^) should not 
be distinguishable (lEj . This argument is supported by 
our numerical results which cannot differentiate the FSS 
scaling in the static and UCM networks. 

Now we move to a typical model exhibiting a nonequi- 
librium phase transition, namely the directed percolation 



(DP) system 0. It has been well known that most of 
nonequilibrium models showing an absorbing-type phase 
transition belong to the DP universality class. Among 
such models, we here consider the contact process (CP) 
and the susceptible- infected-susceptible (SIS) model Q. 

The CP is an interacting particle model on a lattice. 
A particle creates another particle in one of its neigh- 
boring sites with rate p and a particle annihilates with 
rate 1. In the SIS model, the particle creation is at- 
tempted in all neighboring sites. A particle-particle in- 
teraction comes in through disallowance of multiple oc- 
cupancy at a site. As p increases, the system undergoes a 
phase transition at p c from a quiescent vacuum (absorb- 
ing) phase to a noisy many-particle (active) phase in the 
steady state. Near the absorbing phase transition, the 
order parameter (particle density) p ~ e 13 , the fluctua- 
tions x' — N(Ap) 2 ~ e~ 7 , the susceptibility x ~ l e l~ 7 j 
the correlation length £ ~ the relaxation time 

t ~ | e | ~ ^* , and the survival probability P s ~ e 13 with the 
reduced coupling constant e = (p — p c )/p c - It is known 
that (3 = f3' due to the time-reversal symmetry in the 
DP systems (i~9j and 7' 7^ 7 in general nonequilibrium 
systems. 

Consider the droplet (cluster) excitation starting from 
a localized seed in the absorbing phase. The average 
space-time size S of a cluster is estimated as 

S~r e ^~\er, (7) 

where Tg and £ c are the average lifetime and typical size of 
a droplet, respectively. Usually Tg diverges near the tran- 
sition as t £ ~ \e\- v * +f) ' for u t > f3' [H, but n is a 0(1) 
constant otherwise. In the MF regime, it is shown later 
that the latter always applies. The droplet size diverges 
as £ c ~ |e| ~ v t , which leads to a — dv T + maxj^t — /3', 0}. 
It is well known that the susceptibility is proportional to 
the cluster mass, which yields 7 = a — (3 [19| , Finally, we 
arrive at the generalized exponent relation as 

1 = dv T - /3+ max{i/ t -0',Q}. (8) 

The fluctuation exponent 7' satisfies the standard hyper- 
scaling relation as 7' = dv T — 2(3. 

In SF networks, we propose a phenomenological modi- 
fication of the MF Langevin equation describing the DP 
models, similar to the free energy modification of the 
Ising model in Eq. : 

j t p{t) = ep- bp 2 - dp - 1 + y/pqit), (9) 

where p(t) is the particle density at time t and 7/(i) is 
a Gaussian noise. Our modification to the standard MF 
theory comes in by the third p ^ 1 term and it is straight- 
forward to show that the exponent 6 = A for the CP and 
9 = A - 1 for the SIS. 

By dropping the noise term, one may easily get the 
MF steady-state solution for p. We find that j3 = 1 for 
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FIG. 2: (Color online) Data collapse of p for the UCM net- 
work with A = 2.75, using p/v = 0.58 and v — 2.4. In- 
set: double logarithmic plots of the critical decay of p and x' 
against N with slopes f3/u = 0.58(1) and j'/u = -0.16(2). 



networks, where no cutoff dependence is found. The SIS 
model shows bigger finite-size corrections, which hinder 
us to estimate the transition point accurately. Neverthe- 
less, our estimates for the steady-state exponent ratios 
are consistent with our predictions within errors (lij . 

In summary, we have explored the finite-size scaling 
(FSS) behavior in complex networks. Based on the 
droplet-excitation argument, we conjectured the FSS ex- 
ponent values, which have been confirmed by extensive 
numerical simulations on the static and UCM networks. 
Real networks found in many biological, economical, and 
social systems are highly heterogeneous and also quite 
small in size. Our results will provide the essential infor- 
mation on analyzing the data on these networks. 

This work was supported by Korea Research Foun- 
dation Grant funded by the Korean Government 
(MOEHRD) (KRF-2006-331-C00123) and by the BK21 
project. 



6 > 3 and /3 = 1/(9 - 2) for 2 < 9 < 3. For 9 < 
2, there is no phase transition at finite p. The same 
result may be obtained from the well-established the k- 
dependent noiseless MF theory Moreover, the time- 
dependent solution leads to vt = 1 and the inclusion of a 
weak external field leads to 7 — 1 for all 9 > 2. Notice 
that v t <(3= 13'- 

Utilizing the exponent relation of Eq.©, we summa- 
rize the results for j3, 7, and the FSS exponent v = dv T — 
1 +/3 as 



(/3,7>*) 



1 -1 6-1 

0-2> L > 6-2 

1, 1, 2 



for 2 < 9 < 3, 
for 9 > 3, 



(10) 



where 9 = X for the CP and 9 = X — 1 for the SIS. 
For both cases, ut = 1, cr = v, and 7' = 1 — j3. As 
in the Ising model, a naive power counting for the local 
Langevin equation with the spatial fluctuation term V 2 p 
yields the same result for D by using v = d u v G [l6j ]. 

We performed numerical simulations for the CP on the 
UCM networks at various values of A up to N = 10 7 . We 
measure p and x' = A(Ap) 2 . The transition point p c 
is estimated by the power-law temporal dependence of p 
and also by the power-law size dependency of p in the 
steady state. At criticality, the temporal dependence is 
given as p ~ t~PI Vt and, in the steady state, p ~ N~P/ U 
and x' ~ A 7 By collapsing the off critical data in the 
steady state, we can estimate the value of v. 

In Fig. [21 the data collapse is shown for A = 2.75 in 
the UCM networks. Our estimates for the exponents 
summarized in Table I agree well with our predictions, 
in general. Our data become a little bit weaker close to 
A = 2 where high heterogeneity in networks yields big 
corrections to scaling. We also confirmed the relation of 
a = v by investigating the external field effect at crit- 
icality. Simulations were also performed on the static 
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